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I Abstract 

' We suggest an alternative mathematical model for the massless neu- 

trino. Consider an elastic continuum in 3-dimensional Euclidean space 
and assume that points of this continuum can experience no displace- 

1 I ' ments, only rotations. This framework is a special case of the so-called 
O , Cosserat theory of elasticity. Rotations of points of the continuum are 

described by attaching to each point an orthonormal basis which gives 
}_( ' a field of orthonormal bases called the coframe. As the dynamical vari- 

ables (unknowns) of our theory we choose a coframe and a density. We 
write down a potential energy which is conformally invariant and then 
I ^ incorporate time in the standard Newtonian way, by subtracting kinetic 

^ . energy. Finally, we rewrite the resulting nonlinear variational problem in 

QQ ' terms of an unknown spinor field. We look for quasi-stationary solutions, 

, i.e. solutions that harmonically oscillate in time. We prove that in the 

CNj ■ quasi-stationary setting our model is equivalent to a pair of massless Dirac 

' equations. The crucial element of the proof is the observation that our 

, Lagrangian admits a factorisation. 

o : 

O ■ 1 Introduction 

. ^ I The massless Dirac equation is a system of two homogeneous linear complex 

' partial differential equations for two complex unknowns. The unknowns (com- 

^ I ponents of a spinor) are functions of time and the three spatial coordinates. 

■ - - I This equation is the accepted mathematical model for the massless neutrino. 

The geometric interpretation of the massless Dirac equation is rather com- 
plicated. It relies on the use of notions such as 

• spinor, 

• Pauli matrices, 

• covariant derivative (note that formula ([5]) for the covariant derivative of 
a spinor field is quite tricky). 
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There is also a logical problem with the massless Dirac equation in that it 
predicts the existence of four essentially different types of plane wave solu- 
tions which are called left-handed neutrino, right-handed neutrino, left-handed 
antineutrino and right-handed antineutrino. Only left-handed neutrinos and 
right-handed antineutrinos are observed experimentally. 

The purpose of this paper is to formulate an alternative mathematical model 
for the massless neutrino, a model which is geometrically much simpler. The 
advantage of our approach is that it does not require the use of spinors, Pauli 
matrices or covariant differentiation. The only geometric concepts we use are 
those of a 

• metric, 

• differential form, 

• wedge product, 

• exterior derivative. 

Our model also overcomes the logical problem mentioned in the previous para- 
graph in that it predicts the existence of only two essentially different types of 
plane wave solutions. These correspond to clockwise and anticlockwise rotations 
of the coframe. 

The paper has the following structure. In Section [2] we introduce our nota- 
tion and in Section [3] we formulate the massless Dirac equation. In Section [4] we 
formulate our mathematical model. In Section [5] we rewrite our Lagrangian in 
relativistic form; we do this to show that our Lagrangian looks simpler in rela- 
tivistic notation, though we do not pursue the relativistic approach consistently 
in the current paper. In Section [6] we translate our model into the language of 
spinors. In Section[7]we prove Theorem [1] which is the main result of the paper: 
this theorem establishes that in the quasi-stationary case (prescribed oscillation 
in time with frequency tu) our mathematical model is equivalent to a pair of 
massless Dirac equations. The crucial element of the proof of Theorem [1] is the 
observation that our Lagrangian admits a factorisation; this factorisation is the 
subject of Lemma [2j Section [8] deals with plane wave solutions. The concluding 
discussion is contained in Section [9l 

2 Notation and conventions 

We work on a 3-manifold M equipped with prescribed negative (i.e. with 

signature ) metric g^^. We choose negative metric on the 3-manifold M 

in order to facilitate the subsequent extension (see Section [5]) to a Lorentzian 
metric 



of signature H on the 4-manifold M x M. We denote time by x° and local 

coordinates on Af by x", a = 1, 2, 3. 




(1) 
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All constructions presented in the paper are local so we do not make a priori 
assumptions on the geometric structure of {M,g}. 

Our notation follows [2 [11 IS]- In particular, in line with the traditions of 
particle physics, we use Greek letters to denote tensor (holonomic) indices. The 
only difference with references [TJ [U [3] is that, by default, we assume tensor 
indices to run through the values 1, 2, 3 rather than 0, 1, 2, 3. The index will 
be treated separately. 

Details of our spinor notation are given in Appendix A of [T] . We choose the 
zeroth Pauli matrix to be the identity matrix, 

The defining relations for the three remaining Pauli matrices are 

a^y-' + a^^.a"'" ^2g''f'Sa^, a,/3 = 1,2,3, (3) 

^^^^^""' = 0, a =1,2,3. (4) 

We assume that all our Pauli matrices do not depend on time a;*^. Formulae 
©-([I]) mean that we are effectively working on the 4-manifold R x Af equipped 
with Lorentzian metric ([1]). 

By V we denote the covariant derivative on the 3-manifold M with respect 
to the Levi-Civita connection. It acts on a vector field and a spinor field as 
Vau'' := + T'^a'yv'^ and 

VaC do^C + ^ip"' idea" be + T^cja^bc)^' (5) 
respectively, where da ■— d/dx" and 

T/^af ^ I cry } ^^'^''(^^"575 + 9^9aS - dsga^) 

are the Christoffel symbols. We also denote do := d/dx^. 

We identify differential forms with covariant antisymmetric tensors. Given 
a pair of real covariant antisymmetric tensors P and Q of rank r we define 
their dot product as P • Q := -pPai a^Qpi Pr9"^^^ ■ ■ •S"'''^''- We also define 
\\Pf:=P-P. 



3 The Dirac equation 

The following system of two complex linear partial differential equations on 
R X AI for two complex unknowns is known as the massless Dirac equation: 

z(±aVo+a",^V„r = 0. (6) 

Here ^ is a spinor field which plays the role of dynamical variable (unknown 
quantity) and is a function of time x^ G R and local coordinates (a;^, x^, x^) on 
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the S-manifold M. Summation in ^ is carried out over a = 1,2,3. The two 
choices of sign in ([6]) give two versions of the massless Dirac equation which 
differ by time reversal. Thus, we have a pair of massless Dirac equations. 
The corresponding Lagrangian density is 

+ (f a",^v„r - r^\^v„f )] ■ (7) 

Note that the massless Dirac equation and Lagrangian are often called Weyl 
equation and Lagrangian respectively. 

4 Our model 

A coframe z9 is a triplet of real covector fields I?-' e T*M, j = 1,2,3, satisfying 
the constraint 

g = -d^ (g) ~ (g) ~ (g) . (8) 

For the sake of clarity we repeat formula ([8]) giving the tensor indices explicitly: 
ga/3 = — t^qI?^ — "i^Q^?! — '&a'^p- Wc assume our coframe to be right-handed, i.e. 
we assume that det > 0. 

Formula ^ means that the coframe is a field of orthonormal bases. Of 
course, at every point of the 3-manifold M the choice of coframe is not unique: 
there are 3 real degrees of freedom in choosing the coframe and any pair of 
coframes is related by an orthogonal transformation. 

As dynamical variables in our model we choose a coframe "(5 and a positive 
density p. These live on the 3-manifold M and are functions of local coordinates 
(x^, x'^,x^) as well as of time G K. 

At a physical level choosing the coframe as an unknown quantity means that 
we view our 3-manifold M as an elastic continuum and allow every point of 
this continuum to rotate, assuming that rotations of different points are totally 
independent. These rotations are described mathematically by attaching to each 
point a coframe (= orthonormal basis). The approach in which the coframe 
plays the role of a dynamical variable is known as teleparallelism absolute 
parallelism = fernparallelismus). This is a subject promoted by A. Einstein and 
E. Cartan [HEUfi]. 

The idea of rotating points may seem exotic, however it has long been ac- 
cepted in continuum mechanics within the so-called Cosserat theory of elastic- 
ity [7]. The Cosserat theory of elasticity has been in existence since 1909 and 
appears under various names in modern applied mathematics literature such 
as oriented medium, asymmetric elasticity, micropolar elasticity, micromorphic 
elasticity, moment elasticity etc. Cosserat elasticity is closely related to the 
theory of ferromagnetic materials [H] and the theory of liquid crystals [HI [10] ■ 
As to teleparallelism, it is, effectively, a special case of Cosserat elasticity: here 
the assumption is that the elastic continuum experiences no displacements, only 
rotations. It is interesting to note that when Cartan started developing what 
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eventually became modern differential geometry he acknowledged [TP that he 
drew inspiration from the 'beautiful' work of the Cosserat brothers. 
We define the 3-form 

T^'^ := A rfi?i + A dd^ A M^) (9) 

o 

where d denotes the exterior derivative. This 3-form is called axial torsion 
of the teleparallel connection. An explanation of the geometric meaning of the 
latter phrase as well as a detailed exposition of the application of torsion in field 
theory and the history of the subject can be found in J12J. For our purposes the 
3-form ^ is simply a measure of deformations generated by rotations of points 
of the 3-manifold M. 

Note that the 3-form ^ has the remarkable property of conformal covari- 
ance: if we rescale our coframe as 

1?^ e'^d^ (10) 
where /i : M — > R is an arbitrary scalar function, then our metric is scaled as 

9af3 e'^''9af3 (H) 

and our 3-form is scaled as 

without the derivatives of h appearing. 

We take the potential energy of our continuum to be 

P{x") :=- f WT^^'W^pdx^dx^dx^ (13) 

We put a minus sign in the RHS of (jl3p because the metric on our 3-nianifold 
M is negative, see beginning of Section [51 As a result, our potential energy is 
nonnegative, as it should be. 

It is easy to see that the potential energy ([T^ is conformally invariant: it 
does not change if we rescale our coframe as (flUl) and our density as 

p h-> e^'^p. (14) 

This follows from formulae ^ and HT'^'^f = h^a'p^^t^d'"' 9^^9''^ ■ 

Thus, the guiding principle in our choice of potential energy (jl3p is conformal 
invariance. 

We take the kinetic energy of our continuum to be 

if(x°) := / WdW^pdx^dx^dx^ (15) 

JM 

where is the 2-form 

■d ■= hi)^ A do^^ A d^d^ + A d^^^). (16) 
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Unlike ([T^ . we did not put a minus sign in the RHS of P3|) : this is because we 
are now squaring a 2-form rather than a 3-form (the number 2 is even whereas 
the number 3 is odd). As a result, our kinetic energy is nonnegative, as it should 
be. 

It is easy to see that the 2-form is, up to a constant factor, the Hodge 
dual of angular velocity, so (fT5)) is the standard expression for the kinetic energy 
of a homogeneous isotropic Cosserat continuum. One should think of a collection 
of identical infinitesimal rotating solid balls distributed with density p. We think 
in terms of balls rather than ellipsoids because of the isotropy, i.e. we do not 
have preferred axes of rotation. 

We now combine the potential energy and kinetic energy to form 
the action (variational functional) of our dynamic problem: 

S -.^ [ {K{x°)-P{x°))dx° = [ L{d,p)dx°dx^dx^dx^ (17) 

where 

mp):^m\' + \\Tn\')p (18) 

is our Lagrangian density. Note that this Lagrangian density is conformally 
invariant in the Lorentzian sense. The latter means that we rescale time simul- 
taneously with the rescaling of the 3-dimensional coframe. 

Our field equations (Euler-Lagrange equations) are obtained by varying the 
action p7p with respect to the coframe -d and density p. Varying with respect 
to the density p is easy: this gives the field equation -I- ||T''''|p — which 
is equivalent to L{'d,p) — 0. Varying with respect to the coframe i9 is more 
difficult because we have to maintain the metric constraint ([8]); recall that the 
metric is assumed to be prescribed (fixed). 

We do not write down the field equations for the Lagrangian density L('i?, p) 
explicitly. We note only that they are highly nonlinear and do not appear to 
bear any resemblance to the linear Dirac equation ([6|). 



5 Relativistic representation of our Lagrangian 

In this section we work on the 4-dimensional manifold R x M equipped with 
Lorentzian metric ([T]). This manifold is an extension of the original 3-manifold M . 
We use bold type for extended quantities. 
We extend our coframe as 

< - (oi) ' (19) 

e=(j5), J = 1,2,3, (20) 

where the bold tensor index o: runs through the values 0, 1, 2, 3, whereas its 
non-bold counterpart a runs through the values 1, 2, 3. In particular, the 0^ in 
formula (|19p stands for a column of three zeros. 



6 



The extended metric ([T]) is expressed via the extended coframe (|20p as 
g = i9° ® i9° - i?^ «) - ® - (g) i?^ (21) 
(compare with ([5])). The extended axial torsion is 

T^'' := i ( i^" Ad< - 1? V - 1?^ A - A rfi?^) (22) 

=0 

where d denotes the exterior derivative on M x M. Formula can be rewrit- 
ten as 

with and defined by (fTCj) and ([5]) respectively. Squaring (^5)) we get 
llT'^'^lp = lli^'lp + IIT'^^'IP which implies that our Lagrangian density (UHl) can be 
rewritten as 

L(,9,p) = ||T-||V (24) 

The point of the arguments presented in this section was to show that if one 
accepts the relativistic point of view then our Lagrangian density (fT5|) takes the 
especially simple form (I24|) . A consistent pursuit of the relativistic approach 
would require the variation of all four elements of the extended coframe which 
we do not do in the current paper. Instead, we assume that the zeroth element 
of the extended coframe is specified by formula . 



6 Choosing a common language 

In order to compare the two models described in Sections [3] and [J] we need to 
choose a common mathematical language. We choose the language of spinors. 
Namely, we express the coframe and density via a spinor field according to 
formulae 

« = , (25) 

p = s^\detg\, (26) 

(^' + = -s-'e<^o.ab^''<^'dce^ (27) 

€ = s-'C^a^ai^ (28) 

where 

(29) 

Note that throughout this paper we assume that the density p does not vanish. 
This is equivalent to the assumption that the spinor field does not vanish. 

Formulae are effectively a special case of those from Section 5 of 

|13| . the only difference being that now the two spinors in the bispinor are not 
independent but related as rj^^ = cr^^j^^". This is hardly surprising as in the 
current paper we work in a non-relativistic 3-dimensional setting so we do not 
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really need bispinors. We also do not really need to distinguish between dotted 
and undotted spinor indices but we retained this distinction in order to facilitate 
comparison with (13] . 

Formulae establish a one-to-two correspondence between a coframe 

•d and a (positive) density p on the one hand and a nonvanishing spinor field 
on the other. The correspondence is one-to-two because for given i9 and p the 
above formulae define uniquely up to choice of sign. This is in agreement 
with the generally accepted view (see, for example, Section 19 in [H] or Section 
3.5 in [TS]) that the sign of a spinor does not have a physical meaning. 

Remark 1 Formulae h25\) - l2^) look somewhat unnatural in that they are assign 
a special meaning to the element i)^ of our coframe. This can be overcome by 
allowing rigid rotations of the coframe, i.e. linear transformations i— > k'd^ 
where k j,k = 1,2,3, is a constant orthogonal matrix with determinant +1. 
Note that such transformations are totally unrelated to coordinate transforma- 
tions. It is well know that axial torsion is invariant under rigid rotations of the 
coframe, hence our model described in Section^ is invariant under rigid rota- 
tions of the coframe. In particular, it is natural to view coframes which differ 
by a rigid rotation as equivalent. 



7 Quasi-stationary case 

For both models, the traditional one (described in Section [3]) and our model 
(described in Section 2]) , we shall now seek solutions of the form 

C{x°,x\x^, x^) = e-''^'=''C{x\x^, x"^) (30) 

where 7^ is a fixed real number. We shall call solutions of the type (|30p 
quasi- stationary. In seeking quasi-stationary solutions what we are doing is 
separating out the time variable a;°, as is done when reducing, say, the wave 
equation to the Hclmholtz equation. 
Substituting pO|) into ([7]) we get 



V\detg\. (31) 



Substituting ^ into (I25l)-(I2ni) and the latter into ^ we get ■& = Ai?^. 
Hence, Hi^lp — ^uj^ and formula (fTS]) becomes 

L(C)- (^||T-f + ^c.2^p (32) 

where 



P = C^\i,C'V\de^\- (33) 

Note that our 3-dimcnsional metric is negative (see beginning of Section [2]), so 
||jTax||2 ^ Q -j-j^g Lagrangian density (|32p may vanish. In fact, as we shall 
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see from the proof of Theorem [T] in the end of this section, it has to vanish on 
solutions of our field equations. 

In order to compare the Lagrangian densities 
explicit formula for T^^ in terms of the spinor field 
the following 



1| and ([32)1 we need an 
This formula is given in 



Lemma 1 We have 



2i 



where 



a/37 



(34) 
(35) 



is the Hodge dual ofT'^^, £123 :— +1. 



Proof Observe that time does not appear in the formula ([9]) for axial torsion 
(summation is carried out over a = 1,2,3). Hence the result we are proving 
should hold for spinor fields ^ with arbitrary dependence on time and not only 
for quasi-stationary ones. As the expression in the RHS of (|34p is an invariant 
we can also temporarily (for the duration of the proof of Lemma [ij suspend our 
convention (Section that our Pauli matrices do not depend on time. 

In order to simplify calculations we observe that we have freedom in our 
choice of Pauli matrices. It is sufficient to prove formula (pi)) for one particular 
choice of Pauli matrices, hence it is natural to choose Pauli matrices in a way 
that makes calculations as simple as possible. Note that this trick is not new: 
it was, for example, extensively used by A. Dimakis and F. Miiller-Hoissen 

We choose Pauli matrices 



'^aab = ■''jab = « 



lab 



a ^2ab 



where 



"jab 



'lab 
^2ab 



^a,abj 



1\ \ 

1 0, 

i 
-i 
1 
-1 



(36) 



(37) 



WO -V/ 

Let us stress that in the statement of Lemma [T] Pauli matrices are not assumed 
to be related in any way to the coframe i?. We are just choosing the particular 
Pauli matrices p6|) . (|37p to simplify calculations in our proof. 

Substituting §6^, ^ into ([28]) we see that the system (I25|)-(l29l) can 
be easily resolved for ^: solutions are spinors with = and (^-^ which is 
nonzero and real. Thus, we have 



(38) 
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where ft, : M ^ K is a scalar function. Formula ([55)1 may seem strange: we are 
proving Lemma [1] for a general nonvanishing spinor field ^ but ended up with 
formula psp which is very specific. However, there is no contradiction here 
because we chose Pauli matrices specially adapted to the coframe and, hence, 
specially adapted to the corresponding spinor field. 
Formulae (O and ([55)1 imply 

= ^<^"ai'^/3"'(5a^^^ + r^^a^ic) + . . . = ^fT^i^/j'^^V^a'^i, + . . . 

o o 

where s — £,'^cr° ^.h^^ — is the scalar ([25)) and the "dots" denote purely 
imaginary terms. We now write down the spinor summation indices explicitly: 

Finally, substituting explicit formulae ([35]), ([37]) for our Pauli matrices and using 
the "dots" to absorb all purely imaginary terms we get 

= ^ V„1$|2/3 „ ^3a^2 ^ ^1"^2 y^^3/3 

8 

= A . dd^ + A • dd^ + (?92 A d^) ■ d^^] + . . . 

8 

Hence, 

= ^[{d^ A d"^) ■ dd^ + A d^) ■ dd^ + {d^ A d^) ■ d^^]. (39) 
Axial torsion is defined by formula ^ whereas 
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so formula (|35p can be rewritten as 

1 



1 



A z?^) • di}^ + A • di}^ + A z?^) . d,?!]. (40) 
o 

Here we used the fact that because of the negativity of our metric (see (|5])) we 
have f = -1, j ^ 1,2,3. 

Formulae (gO]), ([391) and ([Ml) imply □ 



We are now in a position to establish the relationship between the Lagrangian 
densities ^ and ((32l) . 



Lemma 2 In the quasi- stationary case i30\) our Lagrangian density iSS^) fac 
torises as 

32c. L+JC)L^.JC) 



9 L+,(C)-LDi,(C) 



(41) 



Let us emphasise once again that throughout this paper we assume that the 
density p does not vanish. In view of formulae PT|) . in the quasi-stationary 
case this assumption can be equivalcntly rewritten as 

LUO 7^ iDir(C) (42) 

so the denominator in (j4ip is nonzero. 

Proof of Lemma [2] In the quasi-stationary case (|30p formula takes 
the form 

(*r-)p = I [(f a"^^v„r - ra",i,v„^"i)] v^d^ 



because the factor e 



16 



cancels out. Hence (recall that our metric is negative) 



Substituting this into ((32|) we get 
16 



i(C) = - 



9p 



1 2 



2 2 
CiJ p 



(43) 



Formulae gS]), (I3T|) and 133]) imply gl]). □ 

The following theorem is the main result of our paper. 

Theorem 1 In the quasi- stationary case i30\) a spinor field is a solution of 
the field equations for the Lagrangian density i(C) if and only if this spinor field 
is a solution of the field equations for the Lagrangian density L^[j.{C) or the field 
equations for the Lagrangian density i^j,. (C) ■ 
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Proof Observe that the Dirac Lagrangian densities L^j^. defined by for- 



mula ((3T|) possess the property of scaUng covariance: 

L^:Xe\)^e^''L^.JO (44) 

where /i : M ^ R is an arbitrary scalar function. We claim that the statement 
of the theorem follows from and The proof presented below is an 

abstract one and does not depend on the physical nature of the dynamical 
variable C, the only requirement being that it is an element of a vector space so 
that scaling makes sense. 

Note that formulae (|4T|) and (|44)) imply that the Lagrangian density L pos- 
sesses the property of scaling covariance, so all three of our Lagrangian densities, 
L, Lpjj. and L^^^, have this property. Note also that if C is a solution of the 
field equation for some Lagrangian density C possessing the property of scaling 
covariance then C{C) = 0. Indeed, let us perform a scaling variation of our 
dynamical variable 

C^C + hC (45) 

where /i : M — > R is an arbitrary "small" scalar function with compact support. 
Then = 6 J £(C) = 2 / hC{C) which holds for arbitrary h only if £(() = 0. 

In the remainder of the proof the variations of ^ are arbitrary and not nec- 
essarily of the scaling type (jiSj) . 

Suppose that C is a solution of the field equation for the Lagrangian density 
Lqjj,. [The case when ^ is a solution of the field equation for the Lagrangian 
density i^j^. is handled similarly.] Then L^^^{C,) = and, in view of P^ . 
-^Dir(C) ^ 0. Varying C, we get 

= —J ^L+,^iO = —s J i£ir(C) 

We assumed that C is a solution of the field equation for the Lagrangian density 
-^Dir so (5 / -^Dit.(C) = and formula (06]) implies that 6 J L{C) = 0. As the latter 
is true for an arbitrary variation of C this means that ^ is a solution of the field 
equation for the Lagrangian density L. 

Suppose that C is a solution of the field equation for the Lagrangian den- 
sity L. Then L{() = and formula (|41]l implies that either L^[^.{0 = or 
-^Dir(C) = 0; note that in view of ((42|) we cannot have simultaneously ^^^^(C) = 
and ^011.(0 = 0. Assume for definiteness that L^i^.{0 ~ 0- [The case when 
-^Dir(0 = is handled similarly.] Varying C and repeating the argument from 
the previous paragraph we arrive at (|46)) . We assumed that C is a solution of 
the field equation for the Lagrangian density L so S J L{Q = and formula 
(j46|) implies that S J L'^ii-iO = 0- As the latter is true for an arbitrary variation 



so 
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of C this means that C is a solution of the field equation for the Lagrangian 
density L^^^.. □ 

The proof of Theorem [T] presented above may appear to be non-rigorous but 
it can be easily recast in terms of explicitly written field equations. 

Theorem [T] establishes that in the quasi-stationary case our model reduces 
to a pair of massless Dirac equations ©• There is, however, a small logical flaw 
in this statement. The full time-dependent field equations for our model (as 
well as for the massless Dirac model) are obtained by varying the spinor field ^ 
by a with arbitrary dependence on time whereas in the proof of Theorem [T] 
we have effectively varied the spinor field ^ maintaining quasi-stationarity, i.e. 
we took 

SC{x°,x\x^,x^) ^ e-''^''''SC{x\x'^,x^) (47) 

(compare with ((30|) ). Note that the 6C in the above formula does not depend 
on time. If we now modify (|T7|) so that 5C depends on time this will generate 
an extra term in the field equations, one with the time derivative. It turns out 
that this extra term with the time derivative vanishes. For the sake of brevity 
we do not present the corresponding calculation. 



8 Plane wave solutions 



Suppose that M — M.^ is Euclidean 3-space equipped with Cartesian coordinates 
X = {x^,x'^,x^) and metric gap = diag(— 1, — 1, — 1). Let us choose constant 
Pauli matrices 



'lab 



aab 



'2ab 
■^3ab/ 



I 

'\ 

WO -yj 

(compare with ([571) ) ^'^d seek solutions of the form 



(48) 



(49) 



(compare with (|5D1) ) where w 7^ is a real number, k = (fci,fc2,A:3) is a real 
constant covector and ^ 7^ is a (complex) constant spinor. We shall call 
solutions of the type (|49|) plane wave. In seeking plane wave solutions what we 
are doing is separating out all the variables, namely, the time variable x'^ and 
the spatial variables x = {x^,x'^,x^). 

We look at the field equations of our model described in Section 51 These 
field equations are highly nonlinear so it is not a priori clear that one can seek 
solutions in the form of plane waves. However, plane wave solutions are a special 
case of quasi-stationary solutions and the latter were analysed in Section [71 We 
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know (Theorem [Jl that in the quasi-stationary case our model reduces to a pair 
of Dirac equations Substituting (gH) and ((49]) into (O we get 



±w — fca — fci - ik2\ f C 
-ki + ik2 iw + fca ) I 



0. (50) 



The determinant of the matrix in the LHS of ([SO)) is ui'^ ~ k\ — k"^ — fc| so this 
system has a nontrivial solution C if and only if fc^ + fc| + fc| = tj^ . Our model is 
invariant under rotations of the Cartesian coordinate system (orthogonal trans- 
formations of the coordinate system which preserve orientation), so, without 
loss of generality, we can assume that 



(51) 



Substituting (|5ip into (|50p we conclude that, up to scaling by a nonzero complex 
factor, we have 




C = ( J ) ■ (52) 



Combining formulae ([49]) , ((5T|) and (|52|) we conclude that for each real uo ^0 
our model admits, up to a rotation of the coordinate system and rescaling, two 
plane wave solutions and that these plane wave solutions are given by the explicit 
formula 

^''^ Qe— (-"±-'). (53) 

Let us now rewrite the plane wave solutions (j53p in terms of our original 
dynamical variables, coframe d and density p. Substituting ((21), (1^5]) and (|53p 
into (Elll-lIMl) we get p = 1 and 

\ -sin2uj{x° ±x^)y ^ f cos2w(x" ±x3) I , dl=\o\. (54) 

Note that scaling of the spinor C by a nonzero complex factor is equivalent to 
scaling of the density p by a positive real factor and time shift a;" i— > x° -I- const. 

We will now establish how many different (ones that cannot be continuously 
transformed into one another) plane wave solutions we have. To this end, we 
rewrite formula (1541) in the form 



cos2\uj\{x° + bx'^) \ fasin2\uj\{x^ + bx'^)\ /0\ 

^asm2\uo\{x° + bx^)y cos2\uj\{x° + bx^) j, C = ( Oj (55) 

where a and b can, independently, take values ±1. It may seem that we have a 
total of 4 different plane wave solutions. Recall, however (see Remark [T]), that 
we can perform rigid rotations of the coframe and that we have agreed to view 
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coframes that differ by a rigid rotation as equivalent. Let us perform a rotation 
of the coordinate system 











hi 


X 


xV 




[~X' 



simultaneously with a rotation of the coframe 





/ 






z9V 





It is easy to see that the above transformations turn a solution of the form (|55p 
into a solution of this form again only with 

a —a, b i-^ —b. 

Thus, the numbers a and b on their own do not characterise different plane 
wave solutions. Different plane wave solutions are characterised by the number 
c :— ab which can take two values, +1 and —1. 

The bottom line is that we have two essentially different types of plane wave 
solutions. These can be written, for example, as 

/ cos2\uj\{x° + x^) \ /±sm2\Lu\{x° + x^)\ /0\ 

'dl=iTsin2\u;\{x° + x^)\, -dl^ i cos2\uj\{x° + x^) j , i9l=io\.{56) 

The plane wave solutions describe travelling waves of rotations. Both 
waves travel with the same velocity in the negative x'^-direction. The difference 
between the two solutions is in the direction of rotation of the coframe: if 
we fix x^ and look at the evolution of ((56|l as a function of time then one 
solution describes a clockwise rotation whereas the other solution describes an 
anticlockwise rotation. We identify one of the solutions ((56|) with a left-handed 
neutrino and the other with a right-handed antineutrino. 



9 Discussion 

As explained in SectionlH our mathematical model is a special case of the theory 
of teleparallelism which in turn is a special case of Cosserat elasticity. 

The differences between our mathematical model formulated in Section [1] 
and mathematical models commonly used in teleparallelism are as follows. 

• We assume the metric to be prescribed (fixed) whereas in teleparallelism 
it is traditional to view the metric as a dynamical variable. In other 
words, in works on teleparallelism it is customary to view ([8]) not as a 
constraint but as a definition of the metric and, consequently, to vary the 
coframe without any constraints at all. This is not surprising as most, 
if not all, authors who contributed to teleparallelism came to the subject 
from General Relativity. 
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• We choose a very particular Lagrangian density containing only one 
irreducible piece of torsion (axial) whereas in teleparallelism it is tradi- 
tional to choose a more general Lagrangian containing all three pieces 
(tensor, trace and axial), see formula (26) in [T2|- In choosing our partic- 
ular Lagrangian density ^TE\\ we were guided by the principle of conformal 
invariance. 

The main result of our paper is Theorem [1] which establishes that in the 
quasi-stationary case (prescribed oscillation in time with frequency oj) our math- 
ematical model is equivalent to a pair of massless Dirac equations 

This leaves us with two issues. 

A What can be said about the general case, when the the spinor field ^ is an 
arbitrary function of all spacetime coordinates (x ■^) and is not 

necessarily of the form ([30)1 ? 

B What can be said about the relativistic version of our model described in 
Section El? 

The two issues are, of course, related: both arise because in formulating our 
basic model in Section 3] we adopted the Newtonian approach which specifies 
the time coordinate ("absolute time"). 

We plan to tackle issue A by means of perturbation theory. Namely, assum- 
ing the metric to be flat (as in Section [5]) , we start with a plane wave ([^^ and 
then seek the unknown spinor field ^ in the form 

C{x",x\x^,x^) = e-''^^'="+^-^^C{x'^,x\x\x^) (57) 

where ^ is a slowly varying spinor field. Here "slowly varying" means that 
second derivatives of C can be neglected compared to the first. Our conjecture 
is that the application of a formal perturbation argument will yield a massless 
Dirac equation for the spinor field i^. 

We plan to tackle issue B by means of perturbation theory as well. The 
relativistic version of our model has 3 extra field equations corresponding to 
the 3 extra dynamical degrees of freedom (Lorentz boosts in 3 directions) . Our 
conjecture is that if we take a solution of the nonrelativistic problem which is a 
perturbation of a plane wave (as in the previous paragraph) then, at a perturba- 
tive level, this solution will automatically satisfy the 3 extra field equations. In 
other words, we conjecture that our nonrelativistic model possesses relativistic 
invariance at the perturbative level. 

The detailed analysis of the two issues flagged up above will be the subject 
of a separate paper. 
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